Persistent currents due to point obstacles 
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We discuss properties of the two-dimensional Landau Hamiltonian 
perturbed by a family of identical 5 potentials arranged equidistantly 
along a closed loop. It is demonstrated that for the loop size exceeding 
the effective size of the point obstacles and the cyclotronic radius such 
a system exhibits persistent currents at the bottom of the spectrum. 
We also show that the effect is sensitive to a small disorder. 

1 Introduction 

Persistent currents in rings threaded by a magnetic flux are one of the charac- 
teristic features of mesoscopic systems - see, e.g., [|CGR| , |CWB|| and numerous 



other theoretical and experimental papers where they were discussed. Recall 
that for a charged particle (an electron) confined to a loop C the effect is 
contained in the dependence of the corresponding eigenvalues Ej on the flux 
(j) threading the loop measured in the units of flux quanta, 27rfic|e|"^. Specif- 
ically the derivative dEj/dcp equals —-Ij, where Ij is the persistent current 
in the j-th state. In particular, if the particle motion on the loop is free, we 
have 
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where L is the loop circumference, so the currents depend hnearly on the 
apphed field. 

The above example represents an ideal case when the particle is strictly 
confined to the loop. On the other hand, there are many situations when a 
particle in a magnetic field can be transported being localized essentially in 
the vicinity of a barrier or an interface. A prominent example are the edge 



currents |[Ha| , |MS|| . They attracted a new wave of interest recently when in 
several papers ||BP| , [FGW| , [FM| , |MMP|| conditions were derived under which 
such a transport remained preserved in the presence of a disorder. 

Even weaker interaction giving rise to a magnetic transport, this time 
without a classical analogue, was found in ||EJK]| , where the Landau Hamil- 
tonian was perturbed by a straight equidistant array of 6 potentials. It was 
shown that such an operator has bands of absolutely continuous spectrum 
in the gaps between the Landau levels. An independent evidence for the 



transport in this situation was found in ||Ue| . 

The question we are going to ask in this paper is the following. If we 
take a finite array of identical 6 potentials in the plane making a loop, does 
the mentioned transport give rise to persistent currents? As we will see 
the answer is in general positive if the effective radii of the S potentials 
(specified below - see ( p.4|) ) and the cyclotronic radius corresponding to the 
field intensity are much smaller than the loop size. We also discuss the effect 
of disorder manifested by randomization of the coupling constants and show 
that the effect is sensitive to such perturbations: in an example with a ring 
of 12 obstacles a variation larger than 1% leads to a localization. 



2 The model 

2.1 2D magnetic systems with point obstacles 

As we said above, we are going to study a charged particle in the plane subject 
to a homogeneous magnetic field perpendicular to the plane interacting with 
an equidistant array of point obstacles placed at points Sj, j = 1,. . . ,N, 
belonging to a closed loop. The simplest example is the situation when the 
latter is a circle, 

a, = ("i^cos^, i^sin^") , j = l,...,N. (2.1) 
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For the sake of brevity we shall write a := (oi, . . . , 3^). The obstacles will 
be modeled by 6 potentials so formally the Hamiltonian of the system can 
be written as 
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H = (^-iV - 1(f)) + aj5{x~aj) , (2.2) 

where A is the appropriate vector potential which we choose here in the 
circular gauge, 

A{x) = (--By, -Bx 



To make things simple we have used rationalized units in ( p.2|) putting 
h = c = e = 2m* = 1. On the other hand, the two-dimensional 6 in- 
teraction is an involved object. To give it meaning and to specify the true 
coupling parameter which would replace the formal a/s, we use the definition 
based on self-adjoint extensions which is discussed thoroughly in [|AGHH|| , see 



also | D0 , pHS|| for the situation with a magnetic field. A 6 interaction is 



determined by means of the boundary conditions 

L^{ij,dj)+2najLo{ij,d,) = 0, j = I, . . . , N , (2.3) 
which couple the generalized boundary values 

ipix^ 

Lo{ilj,a) := lim - — — — — - , Li(?/',a):= lim ilj{x) — Lo{ilj,a)\n\x — a\ 
lx-a|^o In |x — a| |x-a|^0L 

We now put a := (ai, . . . , oat) and denote by Hj^^ the Hamiltonian of our 
model, which acts as free outside the interaction support, 
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Kai^i^) = (-*V - 1(f)) ^(f) 



for f ^ a, and its domain consists of all functions which belong to the Sobolev 
class Vr^'^(M^ \ a) and satisfy the boundary conditions ([^.3]). 

The parameters aj characterizing the "strength" of the interaction are 
not trivial |(J!S|| ; this is clear from the fact that the free (Landau) Hamil- 



tonian, which we denote as , corresponds to = oo, j = 1, . . . ,N. The 
difference between aj and the formal coupling constants used in ( p.2|) reflects 
the nontrivial way in which the two-dimensional point interaction arises in 
the limit of scaled potentials. To understand the a/s it is more natural to 
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associate them with the scattering length: by |ES] the point interaction with 
the parameter a reproduces the low-energy behaviour of the scattering upon 
an obstacle of the radius p provided 



a 



2ti 



In p . 



(2.4) 



2.2 Spectral properties of H 



The spectrum of the free operator is the series of infinitely degenerate 
eigenvalues B{2n + 1), n = 0, 1, . . . . A finite number of point interactions 
leaves the essential spectrum intact so aess{H^^) consists again of the Landau 
levels. Moreover, it follows from general principles ||W4 Sec .8.3] that point 
interactions can give rise to at most eigenvalues in each gap of a{H^). It 
is this discrete spectrum which is the object of our interest. 

To find it we can make use of Krein's formula ||AGHH| , App. A] which 
expresses the resolvent kernel of H^g^ by means of that of the free operator, 
specifically 
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;) x') = Gb{x, X; z)+ ^ [^a,a{^)h^'^B{x, Sj] z)GB{a^, x'] z), 

j,m=l 

(2.5) 

where A^ = Jz) is the N x N matrix with the elements 



A 



B 
jm 



[Oj - C{B; z)) 5jm - Gsiaj, a^; z){l-5jm) 



where Gb{x, x'; z) is the resolvent kernel of and ^{B;z) its regularized 
value at the singularity. 



lim 



Gb{x, x'; z) -\ In |x — x'\ 

271 



(2.6) 



which is independent of x. The sought discrete spectrum is given by the 
singularities of the second term at the r.h.s. of ( p.5| ), or in other words, by 
solutions of the spectral condition 



det A^Jz) = . 



(2.7) 



If z = zq satisfies this, it is an eigenvalue of the operator Hj^^ and the 
equation A^g{zQ)d = has a nontrivial solution in C^; the corresponding 
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non-normalized eigenfunction is then 



Af 



B,zo 



(2.8) 



as can be checked by a standard argument |[AGHH|| . 

To make use of the above formulae we need the explicit form of the free 
Green's function. It is known | DMM | to be 



where 



47r 



\B\-z 
2\B\ 



U 



\B\-z 
2\B\ ' 



\B\ 



\x — x 



■l\2 



, (2.9) 



^b{x, x ) := exp 



iB 



\B\ 



"-^{Xi-x\){x2+x'2) - 



and U{u, v; () is the singular confluent hypergeometric (or Kummer) function 
Chap. 13]. Using the asymptotic formula 



T(u) 



[lnC + ^(n)-2^(l) + (0)(C'ln2C)] 



with ip{u) being the digamma function and — '^/'(l) = 7 = 0.577 . . . the Euler 
number we find 



\B\-z 
2\B\ 



27 + In — 



(2.10) 



In this way we have derived the explicit form of the matrix A^^(z) appearing 
in the condition 



3 The results 

In Fig. 1, we show the energy spectra of the system as a function of the 
common coupling strength ai = a. In this example, the Np = 12 of S- 
interactions are placed equidistantly on the ring of circumference L = 2tt. 
We can make this assumption without loss of generality. It is true that 



scaling systems with 2D point interactions reveals peculiar properties ||CS2|| , 
but in effect changing the size of the system is equivalent to the shift of a 
on a constant proportional to logarithm of the scaling factor. One can see 



5 



from this figure that, for our choice of the system size, the couphng becomes 
strong as the parameter approaches the value a = —1, in which case the 
eigenvalues split from the Landau levels of unperturbed system form a bunch 
of negative-energy states; recall that by(p.4|) it corresponds to the scattering 
length p = e~^'^. 

The properties of these states are made clear by plotting their energies 
as the function of the magnetic field B. We set the strength parameter to be 
a = —1. In the top lines of Fig. 2, several lowest states are shown for for three 
cases Np = 6, Np = 12 and Np = 24. As the number of the (5-scatterers 
increases, one clearly sees the emergence of parabolic B dependence, which 
is the hallmark of the persistent current we have been looking for. The 
situation becomes clearer when we plot derivatives dE/dB and d'^E/dB'^ as 
the function of B. A direct evidence of the existence of persistent current 
in our system is shown in Fig. 3, in which the probability flux current is 
depicted for Np =12 case for the magnetic field strengths B = 0.5, B = 1.0, 
B = 1.5 and B = 2.0. 

As with the "usual" persistent currents, one has to take into account the 
imperfection when thinking about the use of the effect in designing various 
devices. Thus the question of the robustness of the persistent current in 
our system is of interest. We test it by allowing random fluctuation of the 
coupling strengths of the 5-scatterers. In Fig. 4, we plot the E, dE/dB and 
d^E/dB"^ as the function of B as before, but letting the coupling strength 
of each scatterer to vary randomly by the amount Aa uniformly around 
the central value a = —1. The fluctuations are taken to be Aa = 0.01, 
Aa = 0.02 and Aa = 0.03 from the left to right. The number of (5-scatterers 
is again set to be Np = 12. Some examples of probability current flow 
patterns at B = 0.5 for each case are shown in Fig. 5. While a more detailed 
analysis of such random perturbations is required. Figs. 4-5 suggest that the 
persistent current generated by the array of 5-scatterers is sensitive to small 
fluctuations of the strength; one typically needs the accuracy of less the 1% 
to have the persistent current intact down to the ground state. 

In conclusion, we have demonstrated a purely quantum mechanism able to 
create persistent currents in systems with arrays of tiny obstacles in magnetic 
field; it is the the wave nature of electrons which is responsible for this 
counterintuitive phenomenon. The result might eventually have a practical 
impact on the design of the quantum ring devices, however, the effect of 
localization by random perturbations deserves a deeper study. 
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Figure 1: The energy spectra of a quantum particle on 2-dimensional plane 
with Np = 12 delta functions of strength parameter a placed equidistantly 
on a ring of radius R = 1 and subjected to perpendicular magnetic field of 
constant strength B = 1. All quantities are in the unit h = c = e = 2m* = 1. 
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Figure 2: The energy spectra of a quantum particle on 2-dimcnsional plane 
with Np delta function of strength parameter a = —1 placed equidistantly 
on a ring of radius R = 1 with perpendicular magnetic field of constant 
strength B. 
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B=1.5 B=2.a 



Figure 3: Profile of probability current for Np — 12 delta-functions of 
strength parameter a placed equidistantly on a ring of radius R = 1 and 
subjected to perpendicular magnetic field of constant strength B — 0.5, 
B = 1.0, B = 1.5 and B = 2.0. 
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Figure 4: The energy spectra of a quantum particle on 2-dimensional plane 
with Np delta function of fluctuating strength parameter a = —1 + Aa 
placed equidistantly on a ring of radius R = 1 with perpendicular magnetic 
field of constant strength B. 
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Figure 5: Profile of probability current for Np — 12 delta- functions of 

fluctuating strength parameter a = —1 + Aa placed equidistantly on a ring 
of radius R = 1 and subjected to perpendicular magnetic field of constant 
strength B = 0.5. 
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